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ABSTRACT 

Imaging a star’s companion at multiple epochs over a short orbital arc provides 
only four of the six coordinates required for a unique orbital solution. Probabil¬ 
ity distributions of possible solutions are commonly generated by Monte Carlo 
(MCMC) analysis, but these are biased by priors and may not probe the full 
parameter space. We suggest alternative methods to characterise possible or¬ 
bits, which compliment the MCMC technique. Firstly the allowed ranges of or¬ 
bital elements are prior-independent, and we provide means to calculate these 
ranges without numerical analyses. Hence several interesting constraints (includ¬ 
ing whether a companion even can be bound, its minimum possible semi-major 
axis and its minimum eccentricity) may be quickly computed using our relations 
as soon as orbital motion is detected. We also suggest an alternative to posterior 
probability distributions as a means to present possible orbital elements, namely 
contour plots of elements as functions of line of sight coordinates. These plots are 
prior-independent, readily show degeneracies between elements and allow read¬ 
ers to extract orbital solutions themselves. This approach is particularly useful 
when there are other constraints on the geometry, for example if a companion’s 
orbit is assumed to be aligned with a disc. As examples we apply our methods 
to several imaged sub-stellar companions including Fomalhaut b, and for the 
latter object we show how different origin hypotheses affect its possible orbital 
solutions. We also examine visual companions of A- and G-type main sequence 
stars in the Washington Double Star Catalogue, and show that > 50 per cent 
must be unbound. 

Key words: astrometry - planets and satellites: fundamental parameters - 
binaries: visual - planets and satellites: individual: Fomalhaut b 


1 INTRODUCTION 

Direct imaging is now a well established technique for 
detecting sub-stellar companions, responsible for the dis¬ 
covery of at least fifty such such objects (exoplanet.eu; 
Schneider et al. 2011). The required contrast sensitiv¬ 


ity favours detection of companions far from their host 
stars, thus complementing the radial velocity and transit 
methods that are more sensitive to smaller orbital sepa¬ 
rations. However an advantage of the latter techniques is 
their ability to measure some orbital elements uniquely. 
Very often such constraints are unavailable for imaged 
companions; their long periods mean observations typi¬ 
cally cover a small fraction of their orbit, so often the 
only known kinematic quantities are their instantaneous 
sky plane position and velocity (e.g. [Biller et al. ;201C)J 
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Neuhauser et al.|2011 |. The line of sight coordinates are 
unknown, so it is impossible to state with certainty that a 
companion is even bound, let alone find a unique orbital 
solution. Instead samples of possible orbital elements are 
generated, typically using a Markov chain Monte Carlo 
(MCMC) analysis to find possible orbital fits to the obser¬ 
vations (e.g. Ford 2006 Chauvin et al.|2012 Kalas et al. 


2013||Beust et ah 2014 Pueyo et al.|2014 |. This method 

has many advantages; it easily incorporates observational 
uncertainties, and can fit orbits using a wide range of ob¬ 
servational constraints (such as additional radial velocity 
measurements, e.g. |Crepp et al. 2012]). However an issue 
arises when considering the resulting samples of possi¬ 
ble orbits, as it is assumed that these samples represent 
the probability distributions associated with the orbital 
elements. Whether they actually do is unclear, as the ar¬ 
bitrary choice of prior probability distributions (hereafter 
“priors”) associated with the unknown quantities influ- 
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ences the outcomes of the MCMC; given the same ob¬ 
servables, different priors will produce different results. 

In this paper we argue that, for companions imaged 
over short orbital arcs, a good alternative to MCMC is to 
show how orbital elements vary with the assumed line of 
sight coordinates. This method has previously been used 
by Golimowski et al. (19981, Golimowski et al. (20001 and 


Hinkley et al. (2010) to characterise the semi-major axis 


and eccentricity of a companion; extending it to all or¬ 
bital elements would allow a reader to extract full orbital 
solutions themselves. It is also prior-independent. We also 
provide general, prior-independent bounds on the orbital 
elements of an imaged companion. Both methods are sim¬ 
ple to use and may be applied as soon as orbital motion is 
detected, and provide complimentary information to the 
MCMC technique. 

The equations and techniques presented in this pa¬ 
per are general, and throughout the paper we apply the 
methods to Fomalhaut b as an example (a companion 
which clearly shows linear motion relative to its host, 
| Kalas et al.|20 13|). The layout of this paper is as follows. 
In Section [2j we outline the problem, and in Section [3] we 
provide a simple method to identify whether a general 
imaged companion can possibly be bound. We discuss 
difficulties in interpreting generated distributions of pos¬ 
sible orbits in Section [4] and how these may be overcome. 
In Section[5]we provide general, prior-independent limits 
on orbital elements. We provide some example applica¬ 
tions of our methods in Section [6] and in Section [7] we 
discuss how the bounds on Fomalhaut b’s orbit change 
when various physically motivated priors are considered. 
In Section [8] we discuss the limitations of our method. 


2 KNOWN AND UNKNOWN 

PARAMETERS 

Consider a binary system, comprising a primary and 
a companion. We assume the total binary mass M is 
known, and we define p = GM where G is the gravi¬ 
tational constant. We assume that multi-epoch imaging 
has been performed over a short orbital arc, yielding three 
more parameters: the binary’s projected angular sepa¬ 
ration Rang, the angular sky plane velocity T4ng of one 
component relative to the other, and the angle tp between 
these projected separation and velocity vectors. Assum¬ 
ing the Earth-binary distance d is known, we can convert 
Rang and I4 ng into projected distance R and velocity V. 
Finally we assume no orbital curvature is observed. See 
Appendix [A] for the derivation of R, V and tp from ob¬ 
servational data. 

We define a coordinate system centred on the pri¬ 
mary with the x, y plane defining the sky plane and 
the projected separation vector lying along the x axis 
(Figure [l]). For simplicity we define tp to lie in the range 
0 < tp < 180°, thus setting whether the observer lies 
at positive or negative z (this will affect angular orbital 
elements, as described in Section [5|. Thus we have two 
unknowns, the companion’s position z and velocity i per¬ 
pendicular to the sky plane, for which we have no priors 
and may only limit by assuming that the binary is bound. 

We define the positive dimensionless parameter B, 



Figure 1 . The coordinate system used in this paper, where 
the (x, y) plane is the sky plane. R and V are the projected 
sky plane position and speed of the companion relative to the 
primary, and tp the angle between these. The z direction is 
chosen such that 0 < tp < 180°, so here points out of the page 
(towards the reader). 


which is a combination of observable (and known) quan¬ 
tities 


B = 


V 2 R. 

2/r 


Equivalently 


( 1 ) 


B = 


f 14. 

87t 2 \pcJ \arcsec/yr 


H 


Ra 


r aag 

arcsec J \^Mq 


M 


( 2 ) 

This parameter will prove useful throughout the paper. 


3 CRITERION FOR A BOUND ORBIT 

Consider a binary with one component at instan¬ 
taneous position r = (R, 0, z) and velocity v = 

{V cos ip, V sin tp, z ) relative to the other. The two objects 
may only be bound if v is below the escape speed, i.e. 
v 2 r < 2/t (where r = |rj and v = |v|). The minimum 
possible values of v and r occur if z = z = 0, hence a 
system may only be bound if its sky plane coordinates 
satisfy 


B < 1. (3) 

Hence Equations [I] and [3] may be used to vet common 
proper motion companions for possible physical associa¬ 
tion (see Section [6]). It follows that for bound orbits 



and 


z 

V 


< Vb - 1 - 1 . 


(4) 


Note that unbound solutions are always allowed, as there 
are no practical upper limits on possible values of \z\ and 


4 DISTRIBUTIONS OF POSSIBLE 
ORBITAL ELEMENTS 

We now consider the companion’s orbital elements. A 
Keplerian orbit is completely described by five elements: 


© 2002 RAS, MNRAS 000, PUTT] 































Orbital constraints on imaged companions 3 


semi-major axis a, eccentricity e, inclination to a refer¬ 
ence plane i, longitude of ascending node (from a refer¬ 
ence direction) and argument of pericentre to. Addition¬ 
ally the true anomaly / determines the body’s location 
along its orbit. We define the x , y (sky) plane as the ref¬ 
erence plane and x (the primary-companion separation 
vector) as the reference direction. Appendix[B]shows how 
these orbital elements are derived from a companion’s 
instantaneous position and velocity. Importantly, all six 
elements require the companion’s six-dimensional coordi¬ 
nates to be known for unique determination. Hence given 
a companion’s sky plane coordinates, we must assume 
values of z and z to generate orbital elements. Varying 
these assumed values generates various orbital solutions, 
given the observed parameters. 

To demonstrate this we require an example compan¬ 
ion which shows linear motion over a short orbital arc, 
and we will use Fomalhaut b. We also use this object as 
our primary example for the remainder of the paper. Fo¬ 
malhaut b is a < Jupiter mass object orbiting a 1.92M 0 


star 7.7 pc from Earth (Kalas et al. 2008; Janson et al. 


2012 Mamajek 2012). Kalas et al. (2013) measured its 


position at four epochs between 2004 and 2012, detecting 
orbital motion but not acceleration. We fit the observed 
positions with a linear trend to derive the sky plane ve¬ 
locity. At the first epoch _R ang = 12.54 ± 0.02 arcsec, and 
we fit Va n g = 0.119 ± 0.006 arcsec/yr and ip = 21 ± 2°. 
Neglecting the observational uncertainties, the resulting 
orbital elements as functions of 0 and i are shown on 
Figure [2] The plot shows that very different orbits are 
possible depending on the unknown line of sight coordi¬ 
nates, and also shows the relationships between orbital 
elements as z and z are varied. We will refer back to 
this plot throughout the paper. Note that whilst obser¬ 
vational uncertainties were omitted for simplicity, they 
need not have been; we could have drawn many combina¬ 
tions of the observables using their measured values and 
uncertainties, and produced a contour plot for each com¬ 
bination. Superimposing these plots would accommodate 
the observational uncertainties by broadening the con¬ 
tour lines, but we choose not to do this in our example 
for clarity. 

Drawing many values of z and z, with uniform priors 
on both variables, results in the orbits shown as the blue 
points on Figure [3] (note only bound orbits are shown). 
Evidently there are many degeneracies between differ¬ 
ent elements that are apparent on these plots. The his¬ 
tograms show the corresponding distributions of possible 
orbital element values, for the assumption of uniform pri¬ 
ors. Again, we have neglected observational uncertainties 
when generating this plot; their inclusion would broaden 
the distributions. 

The histograms on Figure [3] agree with distributions 


of possible orbital elements obtained using MCMC (Kalas 
et al.|[2013 Beust et al.||2014]|. However consideration of 


Figure[2j shows that care must be taken when interpreting 
such histograms. Firstly, if these distributions are inter¬ 
preted as probabilities then it is important to remember 
that these depend on the assumed priors. For example, 
had we chosen a prior which biased towards low \z\ and |z| 
(such as uniform log \z\ and log |z|) then a low inclination 
would be favoured, the opposite of the high inclination 


Uniform z, z Uniform logiofz), logiofz) 


T 



Figure 4. The effect of different priors on the distribution of 
orbital inclination solutions for Fomalhaut b (if bound). The 
left plot shows the distribution found using a uniform prior 
on z and z. The right plot shows uniform logarithmic priors, 
where log 1 g(z/i?) and ]og H} (z/V 7 ) are drawn between 0.01 and 
the maximum values from Equation [4] The dashed lines show 
the maximum allowed inclination, using the method in Section 

m 


inferred from uniform z, z priors. We show this on Fig¬ 
ure [I] While uniform z and z priors may seem the most 
unbiased priors that can be assumed, since it implies no 
knowledge of these variables, this ignorance could be ap¬ 
plied in other ways. For example assuming ignorance of 
the orbital properties, say through a uniform prior on ec¬ 
centricity, will result in a different set of solutions to that 
calculated assuming uniform z and z. Hence the “most 
likely” orbital elements clearly depend on the assumed 
priors. This is further explored in Section [TJ where we 
discuss physically motivated priors that can be applied 
to Fomalhaut b. 

Secondly, distributions of possible orbital elements 
alone say nothing about element degeneracies. For exam¬ 
ple, the semi-major axis and inclination histograms for 
uniform z and z priors on Figure [3] imply that the most 
likely orbit has a ~ 150 au and i ~ 70°. However the a 
versus i plot on Figure [3] shows that no solutions exist 
with both a = 150 au and i = 70°. This is not a problem 
if one has access to the list of orbital solutions used to 
create the histograms, but great care must be taken if 
attempting to interpret possible orbits using histograms 
alone. 

We suggest two methods to characterise possible or¬ 
bits as alternatives to MCMC, which unlike the latter 
technique are both prior-independent and probe the full 
region of allowed parameter space. Firstly, whilst his¬ 
tograms of possible elements are prior-dependent regard¬ 
less of the method used to generate them, the allowed 
ranges of orbital elements are set by the known parame¬ 
ters R, V, p> and y. For example, re-generating Figure [3] 
using different z and z priors produces the same bounds, 
but the density of solutions changes (see Figure]!]). These 
bounds can be calculated (see Section [HJ, so we suggest 
that such bounds are quoted as a prior-independent al¬ 
ternative to histograms of possible solutions. 

Secondly, the degeneracies between elements are also 
prior-independent. Figure [2] fully describes these degen¬ 
eracies, as each combination of z and z corresponds to a 
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Figure 2. The dependence of an imaged companion’s orbital elements on its unknown line of sight position z and velocity z. This 
example is Fomalhaut b, which has -Rang = 12.54 arcsec, \4ng = 0.119 arcsec/yr, <p = 21° and M* = 1.92Mq. Solutions outside 
the dashed lines are unbound; note that a —> oo and e —> 1 as the orbit approaches the unbound limit. The plot does not include 
uncertainties on R, V, p and A7+. which would broaden the contour lines. For Fomalhaut’s observational errors, the la uncertainty 
on the z value corresponding to each contour is of order ±10 au, and that on z is of order ±10~ 2 au yr 1 . Plots like this allow 
the reader to extract possible orbits themselves, by choosing a combination of z and i and reading off the corresponding orbital 
elements from each panel. 


single set of orbital solutions. This figure also allows pos¬ 
sible orbits to be extracted by the reader, without having 
to publish lists of solutions (as would be required with the 
MCMC method). Hence a poorly sampled orbit may be 
defined by contour plots of the orbital elements corre¬ 
sponding to different z and z combinations (as on Figure 
[2|, as well as quoting the aforementioned bounds. These 
two methods would compliment the MCMC method, and 
may even be more informative than the latter (at least 
when motion is only detected over a short portion of the 
orbit, see Section [8]). 


5 ORBITAL ELEMENT BOUNDS 

We now provide general bounds for several orbital ele¬ 
ments, which may be quickly calculated as soon as or¬ 
bital motion is observed without the need for numerical 
analyses. Figures [2] and [3] show that and to may as¬ 
sume any value between 0 and 360°, and this holds for 
any combination of sky plane coordinates. However e has 
a clear minimum, and assuming the companion is bound 
also limits the values of a, i and sometimes /. We provide 
the extreme values of these elements for any bound com¬ 
panion imaged over a short orbital arc (note that defining 
0 < ip < 180° also sets 0 < i < 90° and affects H and lo, 
but makes no difference to a, e and /). 


Firstly examination of Equation |B1| shows that the 
minimum semi-major axis for a bound object will always 
occur if z = i = 0, i.e. 




(5) 


No upper limit may be placed on a, as it tends to infinity 
as the orbit approaches the unbound limit. 

No analytic solution for the minimum eccentricity 
in terms of z and i can be found, so we calculate this 
numerically instead. Firstly we reduce the number of pa¬ 
rameters by making Equations |B1| to |B7| dimensionless 
(i.e. dividing all lengths and velocities by R and V re¬ 
spectively). This yields a natural coupling between the 
variables R, V and p in terms of B; thus we reduce the 
number of parameters from six ( R , V, ip, p, z and z) 
to four ( B , p, z/R and z/V). In Appendix |B| we recast 
the eccentricity equation (Equation B21 in terms of these 
dimensionless variables (Equation |B8| ). The minimum ec¬ 
centricity will occur at a certain combination of z/R and 
z/V, hence its value is only a function of B and ip. On 
Figure [5] we plot the minimum eccentricity contours for 
each combination of B and <p, found by numerically vary¬ 
ing z/R and z/V. This plot is completely general for any 
system imaged over a short orbital arc. 

A numerical fit to the e m i n = 0 contour on Figure 
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Figure 3. Possible orbital solutions for an imaged companion, showing the degeneracies between different elements. Again, 
Fomalhaut b is used as an example. Each point on the off-diagonal plots is a possible orbital solution. Blue solutions were 
generated by drawing z and z uniformly, and the histograms show the distributions of possible orbital elements resulting from this 
prior. For any companion, different priors change the density of solutions within the bounds but not the bounds themselves. Here 
observational uncertainties are omitted; these would blur the plots, but by a small amount in this case. We also show solutions for 
some physically motivated priors specific to Fomalhaut, as detailed in Section [7] yellow points are orbits within 5° of the system’s 
debris disc plane, green orbits have semi-major axes within 10 au of that of the disc, and red orbits pass within 5 au of the disc 
when Fomalhaut b crosses it in projection in ~ 50 years’ time. Note a —> oo as e —> 1, however we only show a < 1100 au for 
clarity. A degeneracy exists between fldisc? ^disc and f^disc + 180°, w^isc + 180° ( |Beust et al.|[2014| ), and we plot both solutions 
here. The dashed magenta lines on the versus cu panels show zu = revise? where vo = Q + lj. 
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Figure 5. The minimum possible eccentricity of a companion 
as a function of its sky plane coordinates. The contours were 
found numerically by varying the line of sight position and 
velocity at each point. This plot is completely general, and 
the point shows B and p for Fomalhaut b as an example. 

[5] shows that a companion cannot have a circular orbit 
unless 

-B < 0.2 + 0.3exp(—0.04|ip — 90°|), (6) 

where p is in degrees. In particular note that the eccen¬ 
tricity cannot be zero if B > 0.5, as for a circular orbit 
v 2 r/p = 1 (recalling that r and v are the relative 3- 
dimensional position and velocity respectively). There is 
no upper limit on eccentricity. 

We now bound the companion’s inclination. By def¬ 
inition this equals 0° if 2 = 2 = 0, and this solution is 
always allowed (unless p = 0° or 180°, in which case the 
inclination must be 90° if the object is bound). The incli¬ 
nation of bound orbits also has a maximum, again set by 
B and p, and we plot these contours on Figure [6] This 
maximum occurs if the companion is just bound, i.e. at 
e « 1. This is because the magnitude of cos* is set by 
h = |r x v| (Equation |B3| - noting that h.z is an observ¬ 
able and thus independent of 2 and 2 ), which is minimal 
at 2 = 2 = 0 and forever increases as \z\ and \z\ increase. 

The true anomaly is sometimes constrained for a 
bound orbit (as it is for Fomalhaut b). Of particular in¬ 
terest is whether an imaged companion could be near 
pericentre or apocentre. We show this on Figure [7] as a 
general function of B and p. If B < 0.2 then a bound 
companion may have any true anomaly; this is because 
e m in ~ 0 (Figure [5]), so all pericentre orientations (and 
hence true anomalies) are energetically allowed. For this 
reason, the region of B and p space where all / values 
are allowed is very similar to the e m in < 0.1 region on 
Figure [5] If B > 0.2 then, depending on p, it may not 
be possible for a bound companion have certain values of 
/. We highlight regions of B, p space for which a bound 



Figure 6. A bound companion’s maximum possible inclina¬ 
tion to the sky plane, as a function of its sky coordinates. This 
plot is completely general. Since we define 0° < p < 180°, i is 
confined to the range 0° < i < 90°. Note that if Equation [3] is 
satisfied, i = 0° is always a possible bound solution if p A 0° 
or 180°. The point shows B and p for Fomalhaut b. 


companion cannot be at pericentre, apocentre or either; 
in the latter case, it is also known whether the compan¬ 
ion is moving towards pericentre or apocentre. This is 
because if the companion is moving away from the star 
in projection (0° < p < 90°), it would be moving from 
pericentre to apocentre if its orbit lay in the sky plane. 
If it is actually moving from apocentre to pericentre, its 
line of sight position / velocity must be reasonably high. 
However these required line of sight coordinates may be 
too large for the orbit to be bound, particularly if B is 
large, in which case it is not possible for the compan¬ 
ion to be moving towards pericentre and be bound. The 
opposite argument is true if 90° < p < 180°. 

Finally although II and w are always unconstrained, 
for many combinations of B and p the companion’s longi¬ 
tude of pericenter zxj = fl + u> is bounded reasonably well. 
For many imaged companions it should thus be possible 
to identify the orientation of pericentre on the sky, which 
may be interesting if other bodies have been detected in 
the system. On Figure [8] we show the value of zu at the 
centre of its allowed range, and also the half width of 
the range. The range of allowed zu values decreases as 
B approaches 1; this is because the ranges of 2 and 2 
resulting in a bound orbit also decreases as B —> 1. For 
an companion with B « 1 to be bound, 2 and 2 must 
be approximately zero; in this case it can be shown that 
zu 2p — 180°. For example it is clear from Figure [8] 
that the projected location of Fomalhaut b’s pericentre 
(if bound) must lie on the opposite side of the star to its 
current projected location, and this will prove important 
in Section [3 
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Figure 7. General constraints on a bound companion’s true 
anomaly /. A bound companion in the / 7 ^ 180° region cannot 
be at apocentre but may be at pericentre, and vice-versa for 
the f yf 0° region. If a companion cannot be at apocentre or 
pericentre, its true anomaly is further constrained inside 0 ° < 
/ < 180° or 180° < / < 360° (note that the actual range of 
allowed / values will lie within these bounds; for example the 
point shows Fomalhaut b, which actually has 25° < / < 133°). 



Figure 8. Allowed values of a general companion’s longitude 
of pericentre ro = f2 + u. For each combination of B and p 
there are a range of possible values of vj. Colours show the 
value of ~ in the middle of the allowed range, i.e. [max(ro) + 
min(to)]/2. Contours show the half width of the range, i.e. 
[max(n) — min(ro)]/2. Much of the region between the 90° 
and 180° contours is dominated by numerical noise and should 
be ignored. The point shows Fomalhaut b, which has 203° < 
ro < 304° (i.e. —157° < m < —124° on the plot). 


APPLICATION TO IMAGED 
COMPANIONS 


Our methods may be applied as soon as a companion’s 
orbital motion is detected, firstly to establish whether it 
could be bound and then to constrain its semi-major axis, 
eccentricity, inclination and true anomaly. Even broad 
constraints are interesting; a non-zero minimum eccen¬ 
tricity may point towards a certain formation mechanism 
or past dynamical evolution (e.g. jMayer et al.|2004|[Chat-| 
terjee et al.|2008|. An unseen mass can make a compan¬ 


ion’s eccentricity appear greater than it really is (Pearce 
|Wyatt fc Ke nnedyp Oldj), hence a high minimum eccen¬ 
tricity could justify a search for other bodies in the sys¬ 
tem. Additionally were a companion’s maximum inclina¬ 
tion incompatible with that derived for another compo¬ 
nent (e.g. a resolved debris disc or the stellar rotation 
axis), one could infer previous scattering or Kozai evo¬ 
lution even if the exact orbit were unknown (e.g. |Ford,| 
Kozinsky fe Rasio|2000 Greaves et al.|2014 |. 


We found the extreme orbital elements of several im¬ 
aged sub-stellar companions as examples, listed in Table 
[I] Each shows linear or near-linear motion relative to its 
host, which is compatible with a bound orbit (B < 1). 
For each example system we drew 10 4 combinations of the 
astrometry, mass and distance measurements, assuming 
Gaussian uncertainties on these values. For each combi¬ 
nation we then derived the corresponding values of R, 
V, ip and B, and also the extreme orbital elements (a m i n 
etc.). Table [l] shows the median values of the resulting 
distributions and the lcr uncertainties. The uncertainties 
on the extreme orbital elements are generally large, as the 
sky motions are typically poorly constrained. Nonetheless 
Fomalhaut b and PZ Tel B (if bound) must be highly ec¬ 
centric. Furthermore a minimum eccentricity compatible 
with zero is still useful for modelling purposes, even if 
the uncertainties are large. Note that we drew 10 4 values 
of the observables to demonstrate the uncertainties on 
the derived quantities; alternatively these errors could be 
roughly estimated by propagating the uncertainties on B 
and p through Equation [5| a nd placing B and ip points 
with error bars on Figures [5] and [6] (as shown for Fomal¬ 
haut b on these figures). 


For a second example we demonstrate how the sim¬ 
ple B < 1 criterion may be used to identify unbound 
companions, by applying the method to visual binaries 
in the Washington Double Star Catalogue (WDS). We 
considered all A- and G-type main sequence (MS) pri¬ 
maries listed in the WDS with parallax measurements 
from Hipparcos (Perryman et al. 19971, with observations 
spanning at least a year; this equates to 1887 visual com¬ 
panions to 1548 A stars and 2352 companions to 2039 G 
stars. We wished to find a lower limit on the fraction of 
these with visual companions listed in the catalogue but 
which simple analysis shows cannot be bound, by calcu¬ 
lating a lower limit on B for each one. From Equation 
[l] B is minimised if R and V are minimised and p is 
maximised. For each companion we found a lower limit 
on V from the difference between its projected position 
relative to the primary at the first and last observational 
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Companion 

d / pc 


/ 

M 0 

R / au 

V / au yr 1 

v/ ° 

B 

min(a) / au 

min(e) 

max(i) / ° 

Ref. 

2M 0103(AB) b 

47 ± 3 

0.36 

± 

0.02 

f S2±| 

tn oc+ 0.08 

u -' 3D -0.07 

46 l? 2 

u - zo -0.03 

55lf 

Q Q + 0.6 
U ' U -0.0 

76l? 0 

1 

Fomalhaut b 

7.70 ± 0.04 

1.92 

± 

0.02 

96.6t°;i 

o qi+ 0.05 
u - yi -0.04 

21+ 2 3 

0 53+ 0 06 
u - o,3 -0.04 

ioill s 

u - dd -0.06 

78l§ 

2 

GJ 504 b 

17.56 ± 0.08 

1.22 

± 

0.08 

43.6t°'J 

f) Q+0- 2 

U,y -0.3 

85±Jg 

o-3ll;? 

2 9 l?e 

0 O”^ 0,4 
u - u -o.o 

74l? 9 

3 

Gliese 229 B 

5.77 ± 0.04 

0.570 

± 

0.002 

44.8 ±0.3 

0.5 ±0.1 

144+1° 

0.3+1;? 

30+2° 

Q Q + 0.5 
U ' U _0.Q 

84+?! 

4 

GQ Lup B 

140 ± 50 

0.7 

± 

0.1 

103±§7 

0.4 ±0.1 

1251] 1 

n q+ 0 .5 

u - u -o.o 

47+f 

Q Q + 0.6 

u - u -o.o 

88l? 7 

5 

HR 8799 b 

39 ± 1 

1.5 

± 

0.3 

68l? 

U ‘ yZ -0.05 

95 ±2 

0 44+ 0,2 
u ’”— 0.04 

571?! 

Q Q + 0.3 

u - u -o.o 

6411 

6 

HR 8799 e 

39 ± 1 

1.5 

± 

0.3 

14.4 ±0.4 

2 Q+0.2 

ZiU _0.3 

9811 

0 4D“*" 0,2 
u -^ u -0.06 

Hi? 

0 f)+ 0,4 
U ' U -0.0 

67l? 5 

6 

PZ Tel B 

52 ±3 

1.13 

± 

0.03 

13.4t° 0 i 

1.7 ± 0.1 


0.3811.H 

111? 

Q K+0.1 

u -°-o.i 

88 ± 1 

5 

TWA 5 B 

44 ± 4 

0.7 

± 

0.1 

87 ±1 

0.52 ±0.07 

1331* 

o-4ll;? 

62l?° 

Q Q + 0.6 
U ' U -0.0 

7E+ 4 

/0 -13 
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Table 1. Sky coordinates of several imaged sub-stellar companions, for which orbital motion but little to no curvature is detected. 
V and ip were derived from a linear fit to the sky positions, and R is the fitted sky separation at the first observation epoch. Note 
inclination is only bounded by assuming the companion is bound. Errors are ler. Observation references: (l) |Delorme et al.|2013| 
^2jjKalns_et L jyJ|201jj_(3}jKuz2jhara etmdJ|20RS]_(4)jGo]jmowski et al.||1998[ (5) |Ginski et al.||2014| plus refs, therein, (6) refs, in 
GozdziewskG^Migaszewski|2(^A (7) Neuhhuser et al!]2(^) T Coordinates relative to AB’s barycentre. The large uncertainties for 


GQ Lup B result from the uncertainty on its distance. The small discrepancies between our PZ Tel B results and those of Ginski 
et al. ( {2014] ) are likely due to the latter allowing for observations of orbital curvature, and also the reason given in Section|8] 


epochs, and took R to be the smaller of the projected 
separati ons at t hese epochs . We used Mamajek’s online 
table^](Pecaut & Mamaj ek|20 131 to infer each primary’s 
mass from its spectral type and, assuming the primary to 
be more massive than any companions, found an upper 
limit on /r by assuming both stars to be the same mass. 
We combine these to give a lower limit on B, and show 
the distribution of these B values on Figure [9] We find 
that of all visual companions to A-type MS stars listed 
in the WDS, at least 60 per cent (1104 objects) cannot 
be bound. The value is 50 per cent (1264 objects) for G- 
type MS primaries. The actual fractions of companions 
that cannot be bound will be higher, as we used con¬ 
servative assumptions for our analysis. A more thorough 
analysis could estimate the companion’s mass from its 
V-band magnitude (assuming it is on the main sequence 
and at the same distance as the primary), or could utilise 
all available astrometry rather than just the first and last 
measurements. However we did not do these here because 
our aim was simply to illustrate the technique using con¬ 
servative assumptions, to show how a simple evaluation 
of B can be used without the need for a detailed orbital 
analysis. 


7 ADDITIONAL Z AND Z CONSTRAINTS: 

AN APPLICATION TO FOMALHAUT 

The value of Figure[2]is also evident when there are other 
physical reasons why 2 and z may be constrained. If re¬ 
gions of 2 , 2 space can be excluded for certain scenarios 
then the resulting allowed ranges of orbital elements can 
be easily read off the contour plots, as well as the degen¬ 
eracies between them. For example, for Fomalhaut there 
is additional information on the geometry of the system 
from images of its debris disc (for this paper we consider 
only the well resolved outer disc, not the warm debris 
component found by Stapclfcldt et al. 2004). Quite how 


1 http://www.pas.rochester.edu/~emamaj ek/EEM_dwarf_ 
UBVIJHK_colors_Teff.txt 



Lower limit on B 


Figure 9. Lower limits on B for companions to all A-type 
(solid line) and G-type (red dashed line) main sequence pri¬ 
maries listed in the WDS, which also have Hipparcos parallax 
measurements. Companions with B > 1 cannot be bound, 
whilst those with B < 1 may or may not be bound depending 
on z and z. 


Fomalhaut b is related to the disc, and also how the com¬ 
panion acquired such a high eccentricity, are the subjects 
of ongoing debate (e.g. |Lawler, Greenstreet fe Gladmanj 


2014 Tamayo|2014 Faramaz et al.|2015 l. In this section 

we consider how some simple hypotheses about Fomal¬ 
haut b’s orbit in relation to this disc affect the allowed 
ranges of its orbital elements. 

Fomalhaut’s outer disc is narrow, hence is well ap¬ 
proximated by a single ring of debris. In our coordinate 
system, a ring tracing the centre of the disc has orbital 
elements aai sc = 141.8 au, edi sc = 0.1, idisc = 66.1°, 


fldisc = 199.2° and eddisc = 30° (Kalas et al. 2013 also 


note the degeneracy between fldisci <Ajisc and fldisc + 180° 
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l.OF 



Figure 10. Values of z and z allowed by some physically 
motivated priors for Fomalhaut b. Yellow points: Fomalhaut 
b orbits within 5° of the disc plane. Green: its semi-major axis 
is within 10 au of the disc’s. Red: its orbit passes within 5 au 
of the middle of the disc when it crosses it in projection. This 
plot may be overlayed on Figure[2]to understand the element 
degeneracies on Figure[3] The degeneracy between fl^iso w^isc 
and 12(ji sc -|-180 o , aj^isc V180° results in two sets of symmetrical 
solutions for each prior. 


Wdisc + 180°). Hence the disc’s longitude of pericenter 
Wdisc equals 229°. For Fomalhaut b we find 203° <vj< 
304°, hence (if bound) its pericentre must be reasonably 
well aligned with that of the disc regardless of its line 
of sight coordinates (as found by Kalas et al. 2013 and 


Beust et al. 2014). This alignment is also clear from the f l 


versus cj panels on Figure [3] where most points he close 
to the zu = tudisc line. 

The most obvious hypothesis is that Fomalhaut b is 
moving in (or close to) the disc plane, as might be ex¬ 
pected for a coplanar planetary system (or one in which 
Fomalhaut b was scattered outwards by a planet copla¬ 
nar with the disc). If Fomalhaut b orbits in the disc plane 
then its present day line of sight coordinates are z = ±73 
au and i = TO. 062 au yr _1 , where the upper or lower 
signs are taken simultaneously and two solutions arise 
from the degeneracy in the disc orientation. Relaxing 
slightly the hypothesis of exact coplanarity so that Fo¬ 
malhaut b’s orbital plane lies within 5° of the disc plane 
(i.e. the angle between the normals of the two orbital 
planes is smaller than 5°) then z and i must lie in the 
yellow regions on Figure [TO] The yellow points on Figure 
[3] show the corresponding orbital solutions, and plot II 
on Figure m shows the resulting distribution of eccen¬ 
tricity solutions. By overlaying Figure |To| onto Figure [2] 
the allowed ranges of (and degeneracies between) the or¬ 
bital elements are completely defined. If Fomalhaut b lies 
within 5° of the disc plane it must have 140 au < a < 
320 au, 0.67 < e < 0.84 and 214° < w < 244°. 

We can consider alternative origins for Fomalhaut b, 
and how these affect its possible orbit. If it originated 
in the imaged debris disc and evolved onto an eccentric 
orbit via secular or resonant perturbations from another 
body, its semi-major axis would be similar to that of the 
disc (prior III on Figure [lTj ). Orbits with semi-major axes 
within 10 au of that of the disc are shown in green on 
Figures [3] and |10[ We can understand these constraints 



0.65 0.75 0.85 0.95 0.65 0.75 0.85 0.95 


e 


e 


Figure 11. The effect of different physical priors on the dis¬ 
tribution of Fomalhaut b’s eccentricity solutions. I: uniform z 
and z priors. II: Fomalhaut b orbits within 5° of the debris 
disc plane. Ill: Fomalhaut b’s semi-major axis is within 10 au 
of the disc’s. IV: Fomalhaut b’s orbit passes within 5 au of 
the disc when crossing it in projection in ~ 50 years’ time. 
The dashed lines show the lower bound on the eccentricity 
calculated using the method in Section [ 5 ] 


on ^ and i using the top left plot of Figure [2j as these 
points follow an a ~ 140 au contour. Again by comparing 
these z and z values to the other plots on Figure [2] we 
see that if Fomalhaut b’s semi-major axis is within 10 au 
of that of the disc then it has 0.7 < e < 0.95 and 17° < 
i < 65°. The small overlap with the yellow points shows 
that it is possible for Fomalhaut b to lie in the disc plane 
and have the same semi-major axis as the disc. Larger 
mutual inclinations are also possible; fl and i define the 
orbital plane, and the fact that many green points lie 
far from the yellow (coplanar) points on Figure |3js n 
versus i panel shows that some orbits have high mutual 
inclinations. The orbital solutions have 203° <za< 234°, 
i.e. an even closer alignment with the disc’s pericentre. 
Pericentre alignment could be expected in this scenario 
because an unseen object perturbing Fomalhaut b would 
also perturb the disc, causing alignment of 
(although this is not always the case: see 
2014). 


these orbits 
Beust et al. 


Finally we consider a third physically-motivated 
prior (prior IV on Figure 111. In ~ 50 years’ time Fo¬ 
malhaut b will cross in projection the ring we use to 
approximate the disc, and we investigate which orbits 
actually pass through the disc when this occurs. This 
scenario could be expected if Fomalhaut b originated in 
the disc, and was later scattered out by another body. In¬ 
terestingly, this scenario may have a directly observable 
consequence, since a physical crossing may increase the 
debris collision rate causing a temporary brightening of 
the region in scattered light (Kal as et al.| 2013). The red 
points on Figures |3| and |To] show orbits passing within 5 
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au of the three-dimensional disc location that Fomalhaut 
b will pass in projection. The linear distribution of red 
points arises because Fomalhaut b’s motion will remain 
roughly linear between its current location and the time 
when it crosses the disc, so every z value has a corre¬ 
sponding z pointing towards the disc (although there are 
two solutions on Figure[lO]due to degeneracy in the disc’s 
line of sight position). For Fomalhaut b to cross the disc 
in ~ 50 years it must be inclined by at least 20° to the sky 
plane, and the allowed ranges of the other elements are 
unconstrained (beyond the general bounds from Section 
[5|. In particular note that the orbit need not be copla- 
nar with the disc to pass through this point; high mutual 
inclinations are also allowed. 

The above discussion shows that specific origin sce¬ 
narios for Fomalhaut b require orbits within restricted 
ranges of the possible parameter space. It also empha¬ 
sised the effect of different priors on the distributions of 
orbital elements. A similar method can be applied to any 
companion if additional constraints can be placed on its 
line of sight coordinates. 


8 LIMITATIONS OF THE METHOD 


The method of characterising an imaged companion’s or¬ 
bit using contour plots and bounds on allowed elements 
has several benefits; it is quick, prior independent, probes 
the full region of parameter space, allows readers to draw 
possible orbits themselves, makes degeneracies between 
elements readily apparent and is particularly useful if 
physical arguments can further constrain z and z (see 
Section [7|. However the method also has limitations, so 
should be used in conjunction with (not instead of) other 
techniques. 

Firstly the method considers only a single position 
and velocity, so it does not account for orbital curvature. 
For a companion with observed curvature, the method 
is still valid at each instantaneous position and veloc¬ 
ity along its path. However many of the resulting or¬ 
bital solutions would be invalid, because the curvature 
would rule them out. That only one position and veloc¬ 
ity is considered is also a limitation even if curvature is 
not observed, as linear motion spanning a long enough 
arc is more constraining than an instantaneous position 
and velocity. An example is PZ Tel B, which has roughly 
linear motion but whose separation has doubled over a 
few years. This extra constraint rules out some of the 
solutions of our simple analysis, hence the possible or¬ 


bital parameters found by Ginski et al. (20141 cover a 


smaller range than our solutions. Finally, incorporating 
additional constraints into our method can sometimes be 
difficult. For example, our method could identify the min¬ 
imum possible eccentricity from astrometry data, but ad¬ 
ditional radial velocity measurements could rule out the 
combination of line of sight coordinates required to give 
this minimum eccentricity. This may not be too large a 
problem though; our method is most applicable to com¬ 
panions far from their host stars, thus radial velocity data 
may not impose significant constraints. 

All of the above limitations are accounted for by 
MCMC, but that method has its own disadvantages and 


does not have all the benefits of our method. There¬ 
fore the two techniques should both be used in tandem. 
In particular our method could be used as soon as or¬ 
bital motion is detected, to quickly calculate the allowed 
bounds on orbital elements. Once curvature is observed, 
an MCMC approach may be more appropriate. 


9 CONCLUSIONS 

Imaging a companion over a short orbital arc yields only 
four of the six coordinates required for a unique orbital 
solution. MCMC analysis is often used to generate prob¬ 
ability distributions associated with the orbital elements, 
but despite the advantages of such methods, the distri¬ 
butions are influenced by the choice of priors, may not 
probe the full parameter space and may be misleading 
if degeneracies are not accounted for. We suggest that 
an alternative way to characterise possible orbits is to 
show how the elements differ as functions of the assumed 
line of sight coordinates, thus allowing a reader to ex¬ 
tract potential orbits themselves. Additionally we pro¬ 
vide a simple criterion a companion must satisfy to be 
bound, as well as general constraints on its semi-major 
axis, eccentricity, inclination, true anomaly and longitude 
of pericentre. Quoting such bounds would be useful when 
characterising possible orbits in a prior-independent man¬ 
ner. Our methods are complimentary to MCMC analysis, 
both techniques having their advantages and disadvan¬ 
tages. In particular our methods may be quickly applied 
on detection of a companion’s orbital motion, without 
the need for numerical analyses. 

As an example we apply our methods to visual com¬ 
panions of A- and G-type main sequence stars in the 
WDS, and find that about 50 per cent cannot be bound. 
We also calculate bounds on the elements of several im¬ 
aged sub-stellar companions, and our results agree with 
literature values found using more complicated analyses. 
Finally we considered the effect of physically motivated 
priors on the allowed orbit of Fomalhaut b; stipulat¬ 
ing that its orbit must lie within 5° of the debris disc 
plane means it must have 140 au < a < 320 au and 
0.67 < e < 0.84. Alternatively if its semi-major axis is 
within 10 au of that of the disc then 0.7 < e < 0.95, and 
if its orbit intercepts the disc when passing it in projec¬ 
tion then i > 20°. 
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of the companion’s Northern N and Eastern E sky off¬ 
sets from the star, we have R = d\JNf + Ef, V = 
£V(N 2 - Nt) 2 + (E 2 - Et) 2 and 


COS ip ■ 


NiN 2 + EiE 2 - N[ - El 

y/N* + El /(iV 2 - Ah) 2 + {E 2 - El ) 2 


• (A2) 


APPENDIX B: ORBITAL ELEMENTS FROM 
CARTESIAN COORDINATES 

The orbital elements of a binary component at position 
r = (x, y, z) and velocity 
other (with h = r x v) are 


sin f2 


■ 1 ix 

hsini 


v = (x, y, z) relative 
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sin / = 


*(l-e 2 L 


he 

where 6 = u> + /, 


r and cos / = 


x(l-e 2 )-' 


(B6) 


r = sgn(r.v) i v 2 — 


h 2 


(B7) 


and the upper signs in Equation |B4| are t aken if h.z > 0 
and the lower if h.z < 0 (Murray & Dermott 19991. Re¬ 
stating the eccentricity equation in terms of the dimen¬ 
sionless parameters B and ip yields 


APPENDIX A: DERIVING R, V AND <p 
FROM ORBITAL DATA 


1 — 4 B (p 2 — 2 pv cos <p + v 2 + sin 2 <p) 
x ((1 + p 2 )- 1/2 - B(1 + v 2 )') 


1 1/2 


, (B8) 


The following equations give R, V and <p in terms where p = z/R and v = z/V. 
of sky separation angle S and position angle P of 
the binary at two epochs (subscripts 1 and 2), sep¬ 
arated by time At. Trivially R = Sid and V = 
yjSi — 2 S 1 S 2 cos(P 2 — Pi) + S$, and <p is given by 


cos <p = 


S 2 cos(P 2 - Pi) - Si 
y/Sl - 2 SiS 2 cos (P 2 - Pi) + Si' 


(Al) 


Note we define 0 < <p < 180°. Alternatively, in terms 
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